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^ . Abstract. We analyze seven year and nine year WMAP temperature maps for signatures of 

5g ; three finite flat topologies ^o = T^, ^i = T^ x M\ and ^2 = S'^xR'^. We use Monte-Carlo 

simulations with the Feldman- Cousins method to obtain confidence intervals for the size of 

the topologies considered. We analyze the V, W, and Q frequency bands along with the 

ILC map and find no significant difference in the results. The 95.5% confidence level lower 

bound on the size of the topology is 1.5Lq for ^q, IALq for ^1, and I.ILq for ^2, where 

Lq is the radius of the last scattering surface. Our results agree very well with the recently 

released results from the Planck temperature data. We show that the likelihood function 

tH- I is not Gaussian in the size, and therefore simulations are important for obtaining accurate 

^D ■ bounds on the size. We then introduce the formalism for including polarization data in the 

analysis. The improvement that we find from WMAP polarization maps is small because of 

the high level of instrumental noise, but our forecast for Planck maps shows a much better 

improvement on the lower bound for L. For the ^q topology we expect an improvement on 

^ . the lower bound of L from 1.7Lq to 1.9Lo at 95.5% confidence level. Using both polarization 

j^ I and temperature data is important because it tests the hypothesis that deviations in the TT 

spectrum at small / originate in the primordial perturbation spectrum. 
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1 Introduction 

The global topology of the universe has been extensively studied in the literature since precise 
measurements of the cosmic microwave background (CMB) radiation became available from 
the COBE and WMAP satellites [1-14]. The latest temperature data from the Planck satellite 
has also been used to constrain the topology of the universe [15]. There are two main methods 
for this study: the circles in the sky method and \ow-l correlations of CMB anisotropics. The 
first method is based on the fact that if the global structure of the universe is smaller than 
the last scattering surface (LSS), the LSS intersects itself resulting in matching circles in 
different directions on the sky. This method has been used to put bounds on the size of 
the universe with the most recent bound of 27.9 Gpc from the WMAP seven year data [13]. 
This method, however, does not work if the global topology of the universe is bigger than 
the diameter of the LSS. The radius of the LSS is Lq = 14.4 Gpc, so the current bound from 
this method is already very close to the maximum value of 28.8 Gpc possibly detectable by 
the method. 

The second main method to detect the topology of the universe uses the correlations 
of the CMB anisotropics on large scales. This method has been discussed in detail in [6], 
where it was also used to analyze three different flat topologies of the universe: ^q = T'^, 
^\ = T^ X M^, and .^2 = S^ x M^. The topology ^i has i infinite dimensions with the other 
3 — z dimensions compactified to the same size L. Ref. [6] found an improvement of about 
20 in — 21n^ for the ^i topology with size 1.9Lo compared to the infinite universe. It is 
therefore of importance to perform other independent tests of these results. 

Previous analyses of the topology of the universe have used only the temperature data 
from CMB experiments, the reason being that the currently available polarization data is 
very noisy. However, polarization data from Planck will have much less instrumental noise 
than currently available data from WMAP, and can be used to improve the constraints 
on the topology of the universe. Signals for non-trivial topology can be interpreted, more 
generally, as evidence that the primordial fiuctuation spectrum deviates from the standard 
prediction of inflationary cosmology at small values of /. A particular scenario, such as the 
torus topology we consider, provides a specific model with a few parameters that can be 
fit to observations. Polarization provides a very important check that anomalies in the TT 



spectrum are primordial, since in this case, deviations in the TT, TE and EE spectrum are 
related as they have a common origin. This would not be the case if deviations in the TT 
spectrum arose from late-time effects after decoupling. 

This paper serves two purposes. Firstly, we develop the formalism that can be used to 
analyze the polarization data in addition to the temperature data to improve the constraints 
on the topology of the universe. Although this may not give a big improvement in the 
results with the currently available polarization data from WMAP, it will be of invaluable 
use with the Planck polarization data that is expected to be released soon. Secondly, we 
use the Feldman-Cousins method [16] to accurately estimate the confidence intervals from 
the likelihood function using simulated sky maps. The previous analysis [6] used the ILC 
temperature map only from seven year WMAP data. We repeat this analysis with a slight 
improvement in the technique using the V, W, and Q frequency bands as well to test if the 
observed effect is based on some anomaly in the ILC map. We also analyze the nine year 
release of these maps to check if there is any improvement in the results from new temperature 
maps. We then analyze the WMAP data using the polarization maps as well, and we give a 
forecast for the Planck data. 

The paper is organized as follows. In Section 2 we describe the likelihood calculation for 
temperature and polarization maps, in Section 3 we discuss the details of the data analysis, 
in Section 4 we describe the usage of the Feldman-Cousins method and how to simulate finite 
sky maps. Our results for WMAP data are summarized in Section 5. We then compare our 
results with the most recent results from Planck temperature data and give a forecast for 
Planck polarization data in Section 6, and we conclude in Section 7. 

2 Likelihood Calculation for Temperature and Polarization Maps 

The likelihood calculation for finite topologies using the temperature data only is described 
in detail in [6]. Here we generalize this analysis to include the polarization data as well. We 
work with the torus topology T^, with sides Li,L2, and L3. The other cases can be obtained 
from this by taking some Lj ^- 00. 

The temperature map T, and polarization maps Q and U are measured by a CMB 
experiment. These maps can be decomposed into spherical harmonics 



T{n)=Y,TimYim{n), (2.1) 

Im 

Q(n) ± iU{h) = Y, T2a/m T2Yim{n) (2.2) 

Im 

where gy/m are the spin-weighted spherical harmonics. The polarization coefficients can be 
further decomposed into real and imaginary parts {E and B modes) [17] 

±2aim = Elm ± iBim (2.3) 

which results in 

Q(n) = IY1 [EU+2Yim{-^) + ^2Yim{n)) + iS,„,(+2llm(n) - _2llm(n))] , (2.4) 

Im 

Uih) = IY1 [EimU2YUn) - .2Yimm + iBimi+2YUn) + .2Yim{n))] . (2.5) 

Im 



All of these maps are, in general, correlated with each other. In order to calculate the 
likelihood .if (m|S')dm, one needs to construct a big covariance matrix (including noise) that 
will include the T, Q, and U modes 

, , , exp r-im*(S' + iV)~iml , ^ 

^H5)dm= ]l;„,,\sJ^,,, ^ ^-n (2.6) 

where m = (T, Q, U), 5 and N are the signal and noise covariance matrices, respectively, Up 
is the number of pixels, and symbol |..| stands for determinant. The signal covariance matrix 
is calculated easier in harmonic space. Let us denote 

{^ImXiim') = ^Iml'm' (2-7) 

where X and Y denote T, E, and B. For an infinite universe these matrices are diagonal in 
harmonic {i, m) space, 

^Iml'm' = ^U'^mm'Ci (2.8) 

while for the finite flat topologies the isotropy of space is broken and these matrices acquire 
non-zero off-diagonal elements [6]. The temperature-temperature correlation matrix has been 
derived in [6]. This result can be easily generalized to include the E and B modes as well 

MfJvm' = {^7r)\-^y^^'J-^^Y.Pd'^^9f{k)9}:*{k)YrJm,^,{k). (2.9) 

k 

where Li, L2, L3 are the sizes of the three sides of the torus topology, P^ is the primordial 
power spectrum of the gauge invariant curvature perturbations C, gf and gf are the radiative 
transfer functions for modes X and Y. The sum over k = {ki, /c2, ^3) runs over the values 

27r , 27r , 27r , , 

h = -f-ni, k2 = -—n2, ks = -—ris (2.10) 

Li L2 -L3 

where ni, 712, n^, are integers. 

All of the symmetry arguments described in [6] for M can be generalized to the other 
modes. In particular, all of these matrices are real, and M-^^, M^^, and M^^ are symmetric. 
From now on we will assume no parity violation, which implies M = 0, M = 0.^ No 
gravitational waves have been detected so far by CMB measurements [18] so we will assume 
that M^^ = as well. 

For low-/ modes, the noise in the temperature map can be ignored, while there is a 
significant noise in polarization maps. For this reason it is easier to calculate the likelihood 
by decomposing it into temperature and polarization parts [17]. This can be done by a change 
of variables. Let us define 

Elm = Elm - ^l'm'lm{^l'm'l"m")~ Ti»m" (2-11) 

which gives 



EimTi1,^>)=0, (2.12) 

\^EimEi,^,j = Mij^ii^i = Mi^ii^i — Mi„^„i^{Mi„j^„i,„^,„) Miiii^iiiii^i . (2.13) 



^CMB lensing introduces non-zero contributions into these matrices, however, lensing efTects are neghgible 
on the large scales we consider. 



We now define 



Q(n) = 2 X] [^«m(+2>1m(n) + -2Ylmin)) + iAm(+2>lm(n) - -2llm(n)) 
Im 

U{n) = lYl [Elm{+2Ylm{il) - ^2Ylm{n)) + iBim{+2Ylm{^) + ^2Ylm{n)) 



(2.14) 



(2.15) 



Then in pixel space we get, (Q{hi)T{hj)) = ,(U{hi)T{hj)) = 0. This ahows for the 
likehhood function to be decomposed into a product of two factors 



^{m\S)dm 



exp 



■^m\Sp + Npy^m 



(27r)"p|5p + iVp|i/2 



exp r-iT*5-iTl 
diii '- ^ ^ — ^dT 

(2^)"p/2|5t|V2 ''■' 



(2.16) 



where m = (Q,U). We have ignored the noise in temperature. The noise matrix for the 
new Q and U variables is the same as for the original Q and U. The new signal covariance 
matrix can be calculated as follows 

Q(nOg(n,))= \ Yl [Mtirn'i+2YUm) + -2YUm)){+2Y,r^,{ii,) + -2Y,r^,in,)) 

Iml'm' 
-M,tf,^,(+2l/m(fli)--2ll„^(A^))(+2l^r^'(flJ)--2llw(A,))] , (2-17) 



Q(n,)C/(n,)) = \ Yl [MB'rn'U2YUil^) + ^2Ylrn{nM+2Yi1rn'{nj)-^2Y{ir.,{n,)) 



Iml'm 
tBB 



Ml'J,^,{+2YUivi)-^2YU■^^)){+2Yll^A■^J) + ~2Ylr„A■^,))\ , (2.18) 



?7(n,)C/(n,)) = -\ Y [Mtf^,{+2YUn,)-.2YUn,)){+2Yir^,{n,)-.2Y,r^,ih,)) 



Iml'm' 
rBB 



-MPJ^,U2Ylm{i^^) + _2l^^„^(A.))(+2l^r^K%•) + _2l^,w(n,))] (2-19) 

where M^^f,^, is given by (2.13). 

Following [17], we further rewrite the polarization part of the likelihood as follows 



^{m\Sp)dm 



exp 



-^{Np^ihyiNp^SpNp^ + NpYHNp^^) 



{27r)^p\Np^SpNp^+Np^\y^ 



\Np\^^dm (2.20) 



which is numerically more tractable since it contains only Np . The temperature part of the 
likelihood calculation is described in detail in [6]. 

3 Data Analysis 

We calculate the radiative transfer functions in (2.9) using the CAMB software [19]. We 
calculate the likelihood function using a modification of the WMAP likelihood code [20-22] 
to include the off-diagonal elements of the covariance matrices in harmonic space. The rest 
of the code is unaffected. Since our analysis is sensitive to large scales only, we use the low-Z 



part of the code which calculates the likelihood in pixel space. The temperature map used is 
the Internal Linear Combination (ILC) map smoothed to 9.183° and degraded to A'sido = 16 
in HEALPIX format [23]. To check the consistency between different frequency bands we 
also do the analysis on V, W, and Q maps separately, smoothed and degraded in the same 
way. The temperature map is masked with the Kp2 mask, after which 2482 pixels are left. 
The noise in the temperature map at that low resolution is negligible, however a 1 jjiK white 
noise is added to each pixel (and a corresponding term to the covariance matrix) to aid the 
numerical regularization of the matrix inversion [24]. As discussed in [25], the noise in the 
low resolution V, W, and Q maps is smaller than 1 ^K meaning that these maps can be 
analyzed in the same way as the ILC map. 

The monopole and dipole terms are removed from the original full sky maps, however 
it is no longer true after applying the mask. These contributions are therefore marginalized 
over by introducing large variance monopole and dipole terms into the covariance matrix. 
It is also necessary to marginalize over the residual foreground contamination. An extra 
parameter ^ is introduced into the likelihood function 

^<^- « = (2.)-,/^(detC)V. "' (-5<^ - «^/)"^-'(^ - «^/') ('•') 

where A is the data vector in pixel space, Aj is a foreground template, C is the covariance 
matrix, Up is the number of pixels. The likelihood is then marginalized over ^ to get 




The polarization map used in the low resolution likelihood code is formed by a weighted 
combination of Ka, Q, and V bands, masked by the P06 mask and degraded to N^\^f, = 8 
[17]. Terms up to /max = 30 are included in the temperature analysis and Zmax = 23 in the 
polarization analysis. It has been shown in [6] that cutting off at / = 30 has no significant 
effect on the temperature analysis and that terms up to / = 20 already contain the essential 
effects of the topology analysis. 

We do a scan over the sizes of the three topologies ^q, ^i, and ^2 with a step of O.lLo- 
These topologies break the isotropy of space, therefore we need to also scan over different 
orientations of the topologies. We describe the orientation of the torus by three Euler angles 
{(f), 9, ip). The notation is the same as in [6]. We denote the coordinate frame of the CMB 
data by x, y, z, and the coordinate frame of the torus by x', y', z' (these axes coincide with 
the axes of the torus. For ^i the z' axis points in the infinite direction, while for ^2 it 
points in the finite direction.). We start by the two frames aligned with each other, then 
rotate the torus counterclockwise around the z-axis by angle (j), then around the new j;-axis 
by angle 9, then around the new z-axis by angle if^. Some of these angles are equivalent 
to each other because of the symmetries of the tori. For each topology and each size we 
need to minimize — 21n^ over the direction. As shown in [6], the likelihood oscillates very 
rapidly as a function of the Euler angles, therefore using a minimization algorithm on the 
whole set of directions may not find the global minimum. For that reason the minimization 
is done in two steps. We first construct a grid in the Euler angle space with step size of 
7r/20 and identify the points in the grid that correspond to equivalent directions of the torus. 
We then scan over all of the non-equivalent directions to find the valleys in the likelihood 
function. By comparing the values of — 21n^ at neighboring points of the grid we check 



that the difference is not greater than 5. We then take the point at which the niininium is 
reached after the initial scan, and all of the points that are above the minimum by not more 
than 5, and do a second scan around these points now using the CERN MINUIT package 
for function minimization [26] . The step size of the first scan is chosen small enough to not 
allow more than one minimum between the neighboring points. This way we can guarantee 
that the second scan will find the true global minimum. 

4 Feldman-Cousins Method and Simulations 

Once the likelihood function is calculated and maximized on the whole parameter space, we 

need to derive the confidence intervals for the size of the topologies analyzed. The likelihood 

depends on the size of the topology in a very complicated way. Namely, the likelihood 

depends on the covariance matrix in pixel space, which depends on the covariance matrix in 

harmonic space, the dependence of which on the size of the topology is given by eq. (2.9). 

The previous analysis [6] used the maximum likelihood method relying on the assumption 

that the likelihood ratios have a gaussian distribution in terms of the size. Although this 

assumption is valid for a large data sample [27], it is not clear if the smoothed low-resolution 

sky map can qualify as a large data sample. In this paper we relax that assumption and 

calculate the confidence intervals exactly by estimating the distribution of the likelihood 

ratios using Monte-Carlo simulations. We use the ordering principle of Feldman and Cousins 

[16] to estimate the confidence intervals. Let us briefly summarize the method in general and 

describe the usage of it for our case. 

In general, the method works as follows. For a fixed point P in the parameter space 

one constructs many simulations Si with the same values of the parameters P. For each 

simulation the best fit point Pbest in the parameter space is found and then the likelihood 

ratio calculated 

^ ^ ^(5.|P) _ ^^_^^ 

■J^[Si\Phcst) 

Equivalently, one can calculate the difference in the logarithms of likelihood 

Ai^ = -2ln^{Si\P) - (-21ni^(5i|Pbcst)) ■ (4.2) 

Then for a confidence level a one calculates a single number A^c such that a of the simulated 
experiments have A^ < A^c- For the real data D the likelihood ratio is calculated for these 
values of the parameters P 

A^D = -2ln^{D\P) - (-21nif(Z)|Pbest)) (4.3) 

and the point P is accepted with confidence level a if A^d < A££c. 

For our analysis we have four new parameters L, (/>, ^, ^, in addition to the standard 
cosmological parameters. The topology analysis is sensitive to the low-/ modes only while the 
standard cosmological parameters are determined from the whole range of /. For this reason 
we fix the standard parameters to their best fit values as found from seven year WMAP data 
[22] and vary the new parameters only.^ We would need to simulate many sky maps for all of 
the different possible values of the four new parameters, however, all of the different directions 
are equivalent for simulations. Namely, we expect the same distribution of likelihood ratios 



The values of the cosmological parameters that we use are lOOfibft — 2.227, Q.ch — 0.1116, JIa = 0.729, 
n, = 0.966, r = 0.085, A?;(0.002 Mpc-^) = 2.42 x 10"^ 



(4.2) for a fixed size L but different directions. For a fixed topology and a fixed size L we 
simulate 500 sky maps for one fixed direction to find A^c for that size and direction, and 
use the same A^c for all of the directions with the same size. 

For the real data we would need to find the likelihood ratio (4.3) for each point of the 
parameter space to decide if that point is accepted or not at a given confidence level. However, 
we are interested in the confidence intervals for the size L only, therefore for each size L we 
maximize the likelhood ratio over the angles and use that value to determine if the size L 
must be accepted or not. This means that a given size L is accepted at a given confidence 
level a if and only if there exists at least one point in the parameter space (L, (p, 9, ip) with 
that value of L accepted at confidence level a. 

The simulations of finite sky maps are performed by first diagonalizing the covariance 
matrices (2.9) in harmonic space, then generating random gaussian samples from the diagonal 
matrices, then switching back to the original basis, and finally Fourier transforming into real 
space using the HEALPIX package [23]. 

5 Results 

We first do the analysis using the seven and nine year WMAP temperature maps. On large 
scales the seven year data already has negligible instrumental noise, so we do not expect a 
significant improvement in the results from nine year temperature data. The best fit points 
for all the different maps are summarized in Tables 1, for the topologies ^q, ^i, and ^2- 
The ^2 topology has one finite dimension (the z axis) and two infinite dimensions, which 
means that it has the continuous symmetry of any rotations around the z axis. For this 
reason, the final rotation by Euler angle '0 does not change the topology, so the likelihood 
function does not depend on the angle. The Euler angles can be converted into the direction 
of the z axis of the topologies in galactic coordinates by the simple prescription, b = ^ — 9, 
and / = </) — ^ . 

In order to obtain confidence intervals for the size L we combine the results from V, W, 
and Q bands, and compare to the ILC results. Since we analyzed sizes with step O.ILq we fit 
a smooth curve to our data points to obtain more accurate results. The rejection confidence 
levels as a function of the size of the topology are shown in Fig. 1 for the seven year and 
nine year WMAP temperature data. The smooth lines are a result of the fit to a 10 degree 
polynomial. For comparison, we show the results from separate maps in Fig. 2 for WMAP 
seven year temperature data for the ^q topology. The results from different maps agree 
reasonably well, however, after combining the V, W, and Q results the agreement with the 
ILC becomes much better. The la and 2a lower bounds for the size of the topologies are 
given in Table 2. The change of results from seven year to nine year data is very small, as 
expected. 

Let us now compare our results with the previous results in [6] where seven year ILC 
temperature map was analyzed. The best fit point for the topology ^i is in close agreement 
with that found previously. For the other two topologies we find points that fit the data 
slightly better than those found previously. Note that the only difference in the analysis up 
to this point is a slightly better scan strategy over the angles. Therefore we find no contra- 
diction to the previously found improvements in likelihood. However, after extracting the 
acceptance confidence limits for the size L of the topologies from the likelihood improvements, 
we obtain confidence intervals that are larger than the ones found before. In particular, the 
infinite universe is accepted at about la confidence level for all three topologies. These 





Figure 1. Plots of exclusion confidence level as a function of L/ Lq obtained from seven year and 
nine year WMAP temperature data only. The upper row corresponds to WMAP7, the lower row 
to WMAP9. The three plots correspond to ^q (left), ^i (middle), and ^2 (right). The green 
triangular points correspond to the ILC map, the blue square points correspond to the combined data 
from V, W, and Q maps. 
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Table 1. Best fit points and improvements in likelihood A^ compared to an infinite universe for the 
topology ^Oj-^i and ^1 from seven year and nine year WMAP ILC temperature data. Note that 
for topology ^1 the likelihood does not depend on the Euler angle ■)/'• 



differences arise because of a completely different method of statistical data analysis used 
in this paper. The previous analysis [6] used the maximum likelihood method which relied 
on the assumption that the likelihood function is nearly gaussian in the size. In Fig. 3 we 
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Figure 2. Plots of exclusion confidence level as a function of LJLq obtained from seven year WMAP 
temperature data only for all of the maps analyzed. The plot corresponds to the ^o topology. 
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Figure 3. Exclusion confidence levels as a function of the improvement in likelihood A^. The 
black solid curve corresponds to the Gaussian approximation, the red dashed curve is found from 
simulations for an infinite universe, the blue dotted curve is from simulations for the .-#o topology 
with size L/ L^ ~ 2.0, the green dash-dotted curve is from simulations for the ^q topology with size 
L/Lo = 1.0. 



plot the exclusion confidence level as a function of the improvement in likelihood (4.2) found 
from the Feldman-Cousins method compared to the theoretical approximation of a Gaussian 
distribution. As we can see, for small sizes {L/Lq = 1.0 for the green dash-dotted curve) the 
Gaussian approximation works better than for larger sizes. For the best fit size of L/Lq = 2.0 
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Table 2. Lower bounds for L/Lq 
WMAP temperature data only. 
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Figure 4. Plots of exclusion confidence level as a function of L/Lo for the ^q topology obtained 
from seven year WMAP temperature data only, temperature and polarization data together, and 
forecast for infinite universes for Planck temperature and polarization data. The green triangular 
points correspond to WMAP7 temperature data only, the blue square points correspond to WMAP7 
temperature and polarization data, and the red round points correspond to the forecast for Planck 
temperature and polarization data. 

and for an infinite universe the curves from simulations are very far from the Gaussian curve. 
This is the reason why simulations are essential in order to find accurate confidence intervals. 
Now we proceed to the results using the polarization data together with the temperature 
data. Since the WMAP polarization data is very noisy we do not expect much improvement 
in the results when we add the polarization data. For this reason we only analyze the 
WMAP seven year ILC map for the ^0 topology only. The best fit point we find using the 
temperature and polarization data together corresponds to L/Lq = 2.0 ± 0.5, (p = 0.446 it 
0.016, 9 = 1.792±0.020, V' = 0.595±0.014. The improvement in -2 In .if is 17.57. Comparing 
to the results for the same map and topology but with temperature data only (Table 1, first 
row) we see that the best fit size stays the same however the direction is different. Also, the 
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Figure 5. Forecasts for exclusion confidence intervals for Planck for the ./#o topology with sizes 1.7Lo 
and 2.0Lo- 

improvement in likelihood for the best fit case compared to the infinite case is less than with 
the temperature data only. 

Fig. 4 compares the exclusion confidence levels as a function of the size of the topology 
for the analysis with and without the polarization data included. As we can see, including 
the polarization data slightly improves the lower bounds compared to the temperature data 
only. Using the polarization data we obtain the lower bound on L/Lq of 1.73 (compared 
to 1.71 from temperature only) at 68.3% confidence level and 1.68 (compared to 1.50 from 
temperature only) at 95.5% confidence level. 

6 Forecast for Planck 

The temperature maps from the Planck satellite have been analyzed to place lower bounds 
on the sizes of the three flat topologies discussed here, among other candidates, in [15]. By 
maximizing the likelihood function over the directions they placed 95% lower bounds on 
L/Lq of 1.66, 1.42, and 1.00 for the ^q, ^i, and .^2 topologies, respectively. Compared 
to our results from Table 2 from the WMAP data we can see that Planck is doing slightly 
better for the ^0 and ^i topologies, but we have a slightly better lower bound for the ^i 
topology than Planck (1.1 instead of 1.0). The analysis methods are similar to the ones used 
in this paper, except the scan strategy over directions. They used 10, 000 randomly chosen 
directions for the likelihood analysis which may not be enough to find the true maximum of 
the likelihood function. As pointed out in [6] the variation of likelihood with the direction is 
very rapid and a fine scan is needed in order to find the true maximum. In the current work 
we used a similar number of directions for our initial scan, but we followed it by a finer scan 
to find the true maximum. 

The main reason why we introduced the formalism for including polarization data in the 
analysis is that the polarization data from the Planck satellite is expected to be released very 
soon, which will contain much less instrumental noise than the WMAP polarization maps. 
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To get an idea of how the Planck data will improve the results we do a simple forecast for 
the ^0 topology. We simulate 100 infinite universes with the noise characteristics of Planck 
[28], find A^if^? from equation (4.3) for all of them, and use the median value to forecast the 
exclusion confidence levels for different sizes. The resulting confidence limits as a function of 
the size are shown in the red dashed curve of Fig. 4. We can see a significant improvement 
in lower bounds that Planck can place compared to WMAP. If the universe is infinite the 
expected lower bound on L/Lq from Planck data is 1.92 at 68.3% confidence (vs. 1.73 from 
WMAP) and 1.89 at 95.5% confidence (vs. 1.68 from WMAP). The lower bounds obtainable 
from Planck polarization and temperature data are also better than the results from Planck 
temperature data alone (1.89 instead of 1.66 at 2a confidence level). 

We also do similar simulations for the ^q topology with sizes 1.7Lq (slightly bigger 
than the lower bound set by Planck temperature data alone in [15]) and 2.0Lo (bigger than 
any lower bound obtained so far from different methods). The results corresponding to these 
cases are shown Fig. 5. Using the temperature and polarization data together the finite 
topology with size 1.7Lq can be detected at 3a level at least (more simulations are needed to 
forecast detection at a higher level) with an error of less than 0.05Lq. The finite topology with 
size 2.OL0 can be detected at about 2a level with an error of 0.3Lo. This proves that Planck 
polarization and temperature data together can put lower bounds on the size of an infinite 
universe comparable to the bounds that have been obtained from the alternative circles in 
the sky test, and can detect finite topologies with sizes greater than those accessible by that 
method. 

7 Discussion and Conclusions 

In this paper we did a detailed analysis of the WMAP seven and nine year temperature 
maps for three flat spatial topologies of the universe. We used Monte-Carlo simulations and 
the Feldman- Cousins method to obtain confldence intervals for the size of the topologies 
analyzed. We further introduced the formalism for including the polarization data in the 
analysis and compared the results with and without the polarization maps for WMAP data 
for the 3-torus topology ^q. We then did a simple forecast for the Planck data. 

One of our main conclusions is that the likelihood function is not Gaussian in the size 
of the topologies considered and simulations are essential for obtaining accurate confidence 
intervals. Although we find very similar improvements in likelihood as a previous analysis of 
the same topologies [6], using simulations changes the confidence levels for the sizes of the 
topologies significantly. 

Using the temperature maps alone from WMAP seven year and nine year data we found 
lower bounds for the size of the topologies analyzed. For the topology .^0 = TT^ the lower 
bound on the size L/Lq at 95.5% confidence level is about 1.5, for the ^i = T^ x M^ topology 
we find a lower bound of about 1.4, and for the ^2 = S^ x M? topology the lower bound 
is about 1.1. The infinite universe is excluded at less than la confidence level for all three 
topologies. We analyzed the V, W, and Q frequency bands in addition to the ILC map 
to make sure that the results are similar. We found no significant improvement from nine 
year WMAP temperature data compared to the seven year data. Our results agree very 
well with the corresponding results from Planck temperature data [15]. The measurement of 
temperature maps on large scales is not improved by Planck compared to WMAP (except 
for better foreground removal), so it is expected to get similar results from these data sets. 
Therefore, our results from temperature maps confirm independently the results from Planck. 
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Including polarization data from WMAP slightly improves the results obtained from 
the temperature data alone. However, a simple forecast shows that the improvement will be 
significant for Planck results (the 95.5% confidence lower bound on the size L/Lq for the ^q 
topology will improve to about 1.9 from 1.7 for an infinite universe, and a finite universe of 
size 2.OL0 can be detected at 95.5% confidence level), so the techniques developed in this 
paper will be very useful for a detailed analysis of the Planck temperature and polarization 
maps. Also, the correlation between the temperature and polarization maps can serve as 
an independent check of the results. If similar signatures are seen in both temperature and 
polarization maps then it will be very unlikely that the effect is a result of some systematics 
in the temperature data or not completely removed foreground. 
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